Enhanced thermopower under time-dependent gate voltage 
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We derive formal expressions of time-dependent energy and heat currents through a nanoscopic 
device using the Keldysh nonequilibrium Green function technique. Numerical results are reported 
for a metal/ dot /metal junction where the dot level energy is abruptly changed by a step-shaped 
voltage pulse. Analytical linear responses are obtained for the time-dependent thermoelectric co- 
efficients. We show that the Seebeck coefficient can be enhanced in the transient regime up to an 
amount (here rising 40%) controlled by both the dot energy and the height of the voltage step. 

PACS numbers: 



Since their discoveries in 1821 by Seebeck [ij and in 
1834 by Peltier [2;], thermoelectric effects have been ex- 
ploited for many applications, such as heat voltage con- 
verters, thermocouples or refrigerators. The Seebeck co- 
efficient, or thermopower S, measures the voltage in- 
duced by a temperature gradient through an open circuit, 
whereas the Peltier coefficient 11 measures the heat flow 
induced by an applied current for no temperature gradi- 
ent. In the linear response regime, the Onsager relation 
gives n = —ST, where T is the average temperature of 
the sample. 

The recent achievement of nanoscale systems has in- 
vigorated research activities in this field and renewed 
the quest of the great thermopower (see Ref. for a 
recent review). On the one hand, stationary Seebeck 
coefficients have been measured in different nanoscale 
systems: quantum dot atomic-size contacts 0, spin 
valves fs], nanowires Q, and carbon nanotube (if- The 
Landauer-Biittiker formalism used for the electrical con- 
ductance was extended to model thermal transport in 
microstructures with many terminals [^, [l3|, including 
inelastic effects [llj. Validity of the conventional ther- 
modynamics linear equations was deeply questioned in 
mesoscopic systems: the Onsager relations between heat 
and charge transport coefficients 0, [l^l , the Wiedemann- 
Franz law which links electrical and thermal conduc- 
tances fl^, and the Fourier law On the other hand, 
time-dependent electric transport also benefits from ac- 
tive research works. Single-electron time-control has 
been demonstrated experimentally 



15, he 



with a fair 

agreement with earlier theoretical developments llTl ISll . 
More recently, deeper issues as memory effec ts 19l| or 
interplay between multiple time-modulations f20l| have 
been addressed theoretically. Besides, calculation of 
time-dependent heat current in a linear phonon chain 
has recently been achieved pl| . However, thermopower 
dynamics lacks both experimental as well as theoretical 
investigations. 

This letter gives a first insight into time-dependent 



nonequilibrium thermoelectric transport. As a major re- 
sult, illustrated in Fig. [1] we show that the thermopower 
can be strongly enhanced during the transient regime in a 
metal/dot /metal device, presented in Fig. [2l Time evo- 
lution of the thermopower exhibits promising features. 
Indeed, it can be significantly modified and controlled by 
changing the dot energy from Eq to eo + 7o ^-t time Iq = 0: 
starting from the stationary value at t < <o, it increases 
during a finite time, and then it converges toward its new 
stationary value at t — > c». 
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FIG. 1: Increase of the Seebeck coefficient in the transient 
regime of a metal/dot /metal junction shown in Fig. (2] £o = 
0.5 with 7o = 0.05 (solid line), 70 =0.1 (dashed line), and 
70 = 0.15 (dotted line). We take symmetric barriers: r_R = 
Fl, to — 0, ef = 0, ksT = 0.1. The unit for energies is F. 

Time-dependent heat current through left (L) or right 
(R) reservoir in equilibrium reads 



(1) 



where {if r/i^)) energy current, (J£ j^{t)) the elec- 

tric current, and fJ-L,R the chemical potential. Time- 
dependent Seebeck coefficient can be obtained from the 
ratio between voltage gradient AV and temperature gra- 
dient AT between the two reservoirs, when both left and 
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right time-dependent electric currents cancel 



Thus, the average energy current reads 



(2) 



whereas time-dependent Peltier coefficient is defined as 
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The general system we consider consists of N energy 
levels in an interacting central region connected to non- 
interacting left and right leads. The total hamiltonian 
reads H 





+ Hr + Hc + Ht with 
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= ^ ^ Vkn{t)cldn + h.C. , 
p—L,R k^p,n 
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where cj, (djj) and Ck (dn) are the creation and annihi- 
lation operators for the leads (dot). Hint is the inter- 
acting part of He (Coulomb interactions, phonons cou- 
pling, etc.). The energy band Ek of the reservoirs can 
be time-dependent through a bias change (source-drain 
voltage). Energy levels of the central region e„ can be 
time-dependent through a modulation of gate voltage. 
For the sake of generality, hopping amplitudes Vkn are 
also allowed to be time-dependent. An extension to mul- 
titerminal system with additional degrees of freedom, e.g. 
spin, is straightforward. In this calculation, we only con- 
sider the electron contribution to the energy current. 
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FIG. 2: (Left) schematic of metal/dot /metal junction with 
current directions for each lead; (right) time-dependence of 
the dot energy level. We define chemical potentials and tem- 
peratures of the left and right leads as: fJ.L,R = ef ± eAV/2 
and Tl,r = T ± AT/2. The Fermi energy is £f and T = 
{Tr + Tl)/2 is the average temperature. 



The energy current operator is related to the time 
derivative of the hamiltonian describing the leads [2^ 
by if = —Hl.r- Calculating commutators [Hl,r,H], 
we end up with (h — I) 



lE,Rit)= i ekit)Vk„{t)cldn + h.c. 

k€L,R,n 
kiEL.R 
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{lE^Rit))= 2Re| Y ek{t)Vkn{t)G<,{t,t)\ 

-Im| Y ^kit)G<,{t,t)\ , (8) 
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where G'^f^,{t,t') — i{c^y{t')ck{t)) is the lead Green func- 
tion, and where we have introduced the mixed Green 
function G'^j.{t,t') = i(c^(t')d„(t)) which obeys the 
Dyson equation 



/oo 

+ G<^,{tM)9l{ti,t') 
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i{d'^,{t')dn{t)) is the dot Green function. 
{t,t') = i/(efc)e-*/t''^*i"'=(*i) and g'^{t,t') = i&{t' - 
are the Green functions of the isolated 
leads. Expression of the energy current becomes 

{lEA^)) = 2Re|Tr{ ^ dh[Q^^{tM)^lR{ti,t) 

^Im|Tr{Ei,;^,(t)G<^(i,0}| , (10) 



where we have defined the self-energy as- 
sociated to energy transfer as H2'^(t,<:') = 

Y.k^L,R^l{t)9l'^{t,t')ek{mk{t'). In these ex- 
pressions, E^^/j and are vectors whereas G^^, ^, 
H^'^ are matrices. Last term in Eq. pUj) is a pure 
reservoir contribution to the energy current. 

We now consider a non- interacting metal/dot/metal 
junction with a single dot level Eq connected to reservoirs 
with constant energy bands. This model is suitable for 
experiments in which Coulomb interaction is weak in the 
dot with strong coupling to reservoirs 111- In that 
case, the energy current takes the form 
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{llAt)) = 2Re| dt 
xTlA^. tut) [G^(t, h)fLAe) + G< (i, h)Q{t - h)] 



(11) 



where O is the Heaviside function, Jl.r. the Fermi- 
Dirac distribution function, pl.r. the density of states, 
and TL,R{ek,t,t') = 2-KpL,R{^k)V^{t)Vk{t') measures the 
strength of the coupling between the dot and each lead. 

We investigate the time-dependent thermoelectric re- 
sponse to a unique change eo(i) = £o + lo{t) with 
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7o(0 = 7o0(^ — ^o)- This models a dot energy switch- 
ing from Eq to £0 + 70 by applying a gate voltage at Iq 
(see Fig. The time-dependent heat current defined 
by Eq. ([T]) for p — L,R is now expressed in terms of the 
spectral function A{e, t) as 
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with [17 1 
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,(13) 



and r = -|- Ffl , for which we assume that Tp does not 
depend on energy, and simply reduces to Fp — 2TTpp\Vp\'^. 




FIG. 3: (Top graph) electric current and (bottom graph) heat 
current through left lead (solid lines) and right lead (dashed 
lines) as a function of time t, for Tr — Tl, to = 0, io ~ 0.5, 
70 — 2.5, ksTL — 1, fcsTi? ~ 0, fJ.L,R = ±0.5. Insets show 
(top inset) the dot occupation, and (bottom inset) the dot 
heat. Dashed lines indicate the stationary limits at t — >■ 00. 

The integration over energy in Eq. ()12p has been per- 
formed numerically. In Fig. |3] is plotted the time evolu- 
tion of electric and heat currents when the dot energy 
is modified abruptly at to ~ 0. Starting from constant 
values at t < 0, left and right currents converge toward 
constant values at t — )■ 00. Between these two limits, 
currents show strong time-dependent variations. Cur- 
rents through the zero-temperature, T/j = 0, right lead (, 
see dashed lines in Fig. [3]) exhibit time oscillations whose 



period is related to £0 and 70, as Eq. explicitly in- 
dicates. These oscillations of the electric current have 
been measured through a Ge dot Concerning the 

heat current, experimental results are still needed. In 
the left lead (see solid lines in Fig. [3]) , these oscillations 
disappear due to thermal effects given by Tl ^ 0. 

Using particle number conservation, the average dot 
occupation number is calculated from electric currents as 
{Nd{t)) ^ e-i Jm,M)dt, where 7J,,(i) = /£(i) + /£(t) 
is the displacement current jl6|. In the top inset of 
Fig- El (Ndit)) globally follows an exponential decrease 
(iVd(oo)}(l - e-*/^-) + {NdiO))e-*/^'- (see dotted line) 
characterized by the relaxation time = h/T: weaker 
is the coupling between the dot and the leads, longer is 
the relaxation time. Time evolution of (iVd(t)) shows os- 
cillations around this decrease, that have been already 
observed in experiments [I61] . 

Eq. Il]) comes from thermodynamic relations in leads 
at equilibrium: dHL n — dQ^j^ + ^J'L,RdNL^ft, where 
Nl,r is the lead occupation number and Ql.r. is the 
lead heat. Similarly, for the dot out-of-equilibrium, 
we write dHc = dQd — ^hdNi, — ^RdNn: the en- 
ergy change in the dot reflects a balance between heat 
variation and charges leaving the dot times their ener- 
gies. Thus, we define and numerically calculate an av- 
erage heat in the dot as {Qd{t)) = J{ll}{t))dt, where 
/^(t) = /jj(i) + I'^{t) - Hrit). This definition perfectly 
agrees with energy conservation including contribution 
from tunneling. In the bottom inset of Fig. [Sj the time 
evolution of A{Qd{t)) = {Qd{t)) — {Qd{0)} shows a be- 
havior similar to the dot occupation number (top inset 
of Fig. [3]). However, dramatic differences occur in the 
stationary regimes, for which {Nd{t)) is always constant: 
the dot heat varies linearly in time, without any violation 
of conservation laws. 

In the linear response limit, the time-dependent See- 
beck coefficient, defined by Eq. can be obtained from 
the approximate Fermi-Dirac distribution function [^: 
/L,ij(e) ~ fo{e) + fl>is)ifXL,R - (e - £f )Ti,fl/T)), where 
/o is the Fermi-Dirac distribution function for the leads 
when hl = I^r- Taking left and right electric currents 
equal zero, we obtain the linear response for S{t) in the 
case of strong coupling to reservoir and small energy vari- 
ation: 



S{t) 



roodef'o{e){e~ep)r{s,t) 
eTC^dsf[,{e)r{e,t) 



(14) 



where T{e,t) = -4FLFfl,Im{A(£, t)}/F is the time- 
dependent transmission coefficient. This result is a gen- 
eralization of the Seebeck coefficient expression obtained 
in the stationary case 0, Q including time-dependence 
of the transmission coefficient. For the steady-state, we 
have (/£) = — (/^) (constant {Nd{t))), as can be seen in 
the top graph of Fig. |31 But in the time-dependent case, 
(/£(t)) = does not imply (/£ (t)) = because of the dis- 
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placement current. Since the Seebeck coefficient is mea- 
sured in an open circuit, we must find the adequate AV 
and AT which simultaneously cancel both currents. It is 
important to emphasize that Eq. (1141) is only valid under 
the following assumptions: i) linear response (small Ay 
and AT in comparison to Eq and T) ; ii) high transmission 
through the barriers (large Tl and in comparison to 
other energies) ; and iii) small gate voltage time- variation 
7o, in comparison to Eq. Indeed, these assumptions allow 
to approximate |yl(£,t)p « —2lin[A{s,t)]/T, and, hence, 
to cancel both {Il{t)) and (/^(i)} at any time. 

Similar expressions can be obtained for the time- 
dependent Peltier coefficient defined by Eq. ^ 



j°-^de{e~eF){fL{e)-fR{s))T{s,t) 



ejr^de{fLie)^Me))Tis,t) 



AT=0 



(15) 



following the same assumptions. In the linear response 
regime, Eqs. (fT4|) and (|15p verify the Onsager relation 
n(<) = -TS{t) at any time. 




FIG. 4: Percentage of increase of the Seebecfc coefficient max- 
imum in the transient regime S'max as a function of 70 /eo 
for ksT = 0.05 (solid line), ksT = 0.15 (dashed line) and 
ksT = 0.5 (dotted line). We take Tr = Fl and £f = 0. 

Fig. [T] is obtained using Eq. (fT4|) of the linear re- 
sponse. It is shown an increase of the thermopower af- 
ter a step-shaped gate-voltage pulse was applied. The 
reason is that in the transient regime, the system is 
much more sensitive to temperature or electrostatic vari- 
ations. Furthermore, we measure the thermoelectric ben- 
efit of the transient regime calculating the percentage 
{S^l^ - Sllf)/Sllf where fi-^^^ is the maximum value 
of S and Sllf = {S{t < to) + S{t oo))/2. In Fig.H we 
have plotted the percentage of thermopower increase as 
a function of 70 /eo- This ratio plays an important role 
since it does control the thermopower increase. In such 
a junction, the transient thermopower can be tuned by 
both io and 70, which depend on the dot structural prop- 
erties and of the applied gate voltage (see Fig. [2]) . Higher 
is the ratio, higher is the thermopower increase. Here, an 
increase up to 40% is obtained at small temperature. 

We have proposed a first approach to heat dy- 
namics in nanoscale junctions. General formula for 



the time-dependent heat and energy currents flowing 
through an interacting resonant-tunneling system have 
been obtained. We show that an enhanced thermopower 
can be generated during the transient regime in a 
metal/ dot /metal junction, and that its maximum value 
can be tuned by both the dot energy and the gate voltage. 
With such numerical investigations, it will be possible to 
go beyond the linear response for the Seebeck and Peltier 
coefficients, and further determine non-linear thermody- 
namic laws. Moreover, we shall consider interacting sys - 
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tems in order to analyze phonon bath contribution 
impact of electron-phonon interaction [ll|, and to study 
influence of charging effects in other experiments [23 |. 
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